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Abstract 

Analytic continuation of viscous shock solution for the generalized Burgers equation with 
polynomial nonlinear source term is investigated. We show that a pertubated wave recovers 
its analyticity in the space variable in the strip limited by the first pair of singularities of the 
wave. 

1 Introduction 

Nonlinear parabolic evolution partial differential equations arise in various fields such as fluids 
mechanics, traffic flow modelling, mathematical biology and economics. Burgers' equation is one 
of the most popular of such type of equation, it was first introduced by Bateman [T] and deeply 
studied by Burgers [4] as a one-dimensional model for turbulence. It received since a constant 
mathematical interest from which Hopf's study [8j may particularly be outlined. Consider the 
Cauchy problem for generalized Burgers' equation : 

df{x,t) d^f{x,t)) _ d^f{x,t) 

dt dx dx^ (1) 

f{x,t = 0) = fa{x) ^ a±,a; -> a±,a+ > a_ 

If $ is a function which satisfies the Gelfand-Oleinik (cf . [5] and [T3] ) entropy condition : 

$(u)-$(q-) 

Vue[a-,a+\, c= < [2) 

a+ — a- a+ — a_ 

this equation enables propagation of Shockwaves f{x — ct + d), solving the ordinary differential 
equation : 

-cf{x) + ^f{x)) + k = iyf'{x), /(±cx3)-a± (3) 

Example of such $ are given by concave function on the interval «_)_]. If in addition, $ satisfies 
the Lax [12] condition, 

$'(«+)< c< (4) 
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the wave is said to be non characteristic, and in this case, the large time asymptotics of the Cauchy 
problem ^ are the shifted waves f{x — ct + d). In our study, $ is assumed to be polynomial. Its 
degree will be denoted n. In the particular case = which will be refered further to as 

the classical case, the wave has the explicit expression : 

fix)=c+ tanh Ux + d)^j—j (5) 

the celerity c and the constant k being constrained by the limits f{±oo) — a± : 

c = 

a_<i>(a+) - a+<i>(a_) 
k = 

— a- 

The shift d may be assumed to be zero up to a space translation of the initial data /o- The aim of 
this note is to figure out how the width of the analyticity strip of the solution /(x, t) is matching the 
width yo{— ^^^^ for the classical case) of the analyticity strip of its asymptotic regime f{x — ct). 
For Navier-Stokes equations, the analyticity radius of the solutions is known to increase provided 
that the initial data lies in a functional space that forces it to be zero at infinity. More precisely, 

• Foias and Temam [5] considered space-periodic solutions and proved that the width of the 
analyticity strip increases like ^/t provided the initial data belongs to the Sobolev space 

and that the H^-norm of the solution remains uniformly bounded in time. They used 
Gevrey-norm estimates. 

• Grujic and Kukavica [lOj proved the analyticity of local-in-time solutions with initial data 
lying in a LP(R")-space and also estimated the growth of the width of the analyticity strip 
as -s/t. 

• Lemarie-Rieusset [T3] has also proved a \/i-growth for solutions with initial data in S"(M^) 
using techniques from harmonic analysis. 

Classical Burgers equation may be seen as a ID version of Navier-Stokes equation. We stress 
that the width of the analyticity strip does not increase indefinitely if the overfall — of the 
Shockwave does not vanish. 

The author wishes to thank G. M. Henkin, his thesis supervisor, to have shared with him this 
quite interesting problem. 

2 Statement of the result 

Theorem 1. Let f{x,t) be a solution of Burgers' equation with initial condition fo{x) such that 
1/0(2:) — /(a:)| —x~>±oo 0(e~"l^l) for some a > and such that J^^ifoi^) ~ f{x))dx = 0. For 
all £ €]0, 1], there exists T* > and M > such that f{x,t) extends holomorphically in variable 

z — X + iy in the strip \y\ < min (^yo{l — s), M^v{t — T*)^ for t > T* , where iiya are the closest 
singularity of the real axis of the wave f 
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Remark 1 : This theorem means that after a transient regime, the width of the analyticity strip of 
a perturbation of a viscous Burgers Shockwave is bounded from below by a non decreasing function 
of time converging monotonically to the width of the analycity strip of the wave, hmited by the 
closest wave's singularity, which is both an infinity and a branching point of order n — 1. This 
singularity is therefore a pole only when n = 2, in the case of the classical Burgers' equation. 

Remark 2 : The proof uses Kato's contraction argument for solving Navier-Stokes equation [TT] . 
some elements of Sattinger's stability theory of waves |15] and estimates on the convergence rate 
in time of the real solution to its asymptotic wave [7] , [H] . 

Remark 3 : At the limit 1^ = 0, the solutions develop a discontinuity called the inviscid Shockwave, 
this phenomenon can be studied as a result of the loss of analyticity, as proposed by Bessis and 
Fournier in [3]- Poles are staying away from the real axis so far the viscosity does not vanish. 



3 Analytic properties of waves 

The celerity c may assumed to be zero, as ^(u) may be replaced by $(u) — cu without loss of 
generality. The viscosity is brought to be equal to the unity after an appropriate scaling of the 
solution. Traveling waves are solutions of an ordinary differential equation : 

/' = <i>(/)-$(a±) (7) 

The analyticity of the solution will now be investigated. 

Proposition 1. Travelling waves solutions f{x) of the Burgers equation are analytic functions of 
z — X + iy except in a set of isolated points where they admit a singularity which is both an infinity 
and a branching point of order n— 1. Moreover, f is analytic on a uniform strip {x + iy : \y\ < yo}, 
where yo is the imaginary part of the closest singularity of f of the real axis. 

Proof. Singularities of solutions of such equations were studied for the first time by Briot and 
Bouquet [5]. A more modern approach may be found in [B]. In the neighbourhood of a regular 
point {zq, fo) G C X C, the problem : 

/' = <i>(/)-$(a±) 

f{zo) = fo 

is known to admit powers series solutions around (^Oj/o) with non vanishing analyticity radius. 
Consider now a point where $(/) = oo, ie / = oo. Let u — j, which yields : 

^=-u2(<f>(^-i)-<l>(a±)) 

az 

or, 

^ ^ I (o-s 

du -u2($(u-i)_$(a±)) ^ 
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Where, 



CU" ^(cq + ClU + 



-w2($(u-i) - $(a±)) 

is an analytic function of u in a neighborhood of u = 0. Hence (u = 0, zq) is a regular point of 
the above equation, which therefore admits an analytic solution z(u) in a neighbourhood of u = 0. 
Differentiating ([8]) w.r.t u gives : 

z 

— =0 , Vl<fc<n-1 
du" \u=o 

And consequently : 

z{u) - zo = """^(70 + 71" + • ■ •) 

And finally : 

f{z) = — - = (z - zo)"^ {So + Si{z ^ zo) + . . .) 

Which is well the expected type of singularity. 
Integrating ([7|), we get that : 

f°° dF 
'° = io WTH^) 
The value of this integral depends on the contour from to infinity : 

zo e Zn + 2i7rResidue ( = — | Z (10) 

where 1 < k < n. Singularities are regularly spaced on lines parallel to the imaginary axis and 
there exists at maximum n = deg($) of those lines. This achieves the proof : yo, the imaginary 
part of the closest singularity of the real axis is positive. 

□ 



Define 



y^iw) = exp (^-i J^' ^'im)dr?j = exp ( -1 £ C^dij 



1 1 



(11) 



\rm-2 mf{ty))-^a+)\i 



Consequently the quantities \f'{iy)\^^^w{iy), \^{f )\^^^w{iy) and are bounded func- 

tion of y in the interval \y\ < yo. 

4 Analyticity of solutions 

4.1 Weighted norm approach 

The function f{x,t) = f{x) + h{x,t) solves (P) if h{x,t) solves : 

dh{x,t) ^ d^'if{x))h{x,t) d'h{x,t) _ dR{f,h) 

dt dx dx"^ dx 
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with initial data ho{x) = f{x) - f{x). R{j,h) = $(/ + h) - $(/) - $'(/). /i. The potential 
H(x, t) — h{rj, t)drj solves the partial differential equation : 

^ + *'(/(.))^ - = -ml f ) (13) 

with initial data Ho{x) = ho{ri)drj, yielding the integral equation : 

H{x,t)^Fo{x,t)- dr G(x,t-T-OR{L^)de. (14) 

where Fq = J^^G{x,t; ^)Ho{^)d^ and G'(x,t;T) is the Gr een function of the linearized operator 
L. Differentiating wrt x, we get the integral equation solved by h{x,t). 

pt p+ca 

hix,t)=d,Foix,t)^ dT Gix,t-T;ORifi^),h{^,T))d^ (15) 

Jq J-oo 

Following Sattinger [TS], the Green function G of the operator L in the above integral equation 
may be expressed as : 



Proposition 2. The Green function G{x,t;^) of the linearized operator L has the expression 

w{x) 



Gix,t;0 = ^K{x,t;0 (16) 



where w stands for the weight function 

w{x) ^ exp (^-^ <f\f{v))dv^ (17) 



and K{x, t; ^) for the Green function of the operator : 

32 



Muix,t) = - - i (^f{x)^"{f{x)) + u{x,t) (18) 

Proof. Let H{x,t) = w~^{x)u{x,t). A direct calculation yields : 
d.,H{x,t)=w-\x) (^^^'if)u{x,t)+dMx,t) 

d^^H{x,t) = W-\x) (^Q<f'(/)' + ^f'{xW{f?j u{x,t)+^'{f)d^:,U + d^,U 

Combining this two relations, we get that LH = w^^MwH. The function H{x, t) = w~^{x)u{x, t) 
is therefore a solution of LH{x, t) — g{x, t) if and only if u{x^ t) solves the linear equation A4{wH) = 
wg. This completes the proof. □ 

Definition. Denote m{y,t) — sgn{y) min{\y\, AI\/t) and Ac = {x + iy : \y\ < c}. Introduce the 
following functional spaces : 
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• the weighted Hardy Space 

K,a,bi^c) = {u€ 0(Ae) : u{R) C R, Mnl^ ^^A^^ = 

sup J^^M^f r \vj(^ + iy)\.\u{^ + iy,t)\^dcY' <Oo\ 
0<y<c (^c- |y|) Vi-oo / J 

wzi/j norm IMIhj^^ ja,)- 

• the Banach space Be,T of holomorphic functions f{.,t),0 < t < T in the strip ^rn(yo(i-e),t)) 
with norm \\-\\e,T defined as 

where, 

2 5 

On = 7 

n — 1 

The forthcomming proposition is also inspired by Sattinger with some modifications: 

• Estimates of the analytic continuation of the Green function K in spaces variables are needed. 

• These estimates are realized in a weighted L^-Hardy norm as introduced in the previous 
definition and not in a first order weighted Sobolev in supremum norm as in Sattinger's 
work. 

Proposition 3. The kernel K{x, t; ^) is analytic in variable x and ^ in the strip A^^^^^ Ms/t) ^''^'^ 
satisfies the following estimates : 

K{x + iy, t,^ + ivi)w{i + ivi)h{£, + irifd£,\\L-^i^ax) 



(20) 



^ (a-ri)^ const.e 1 , ,,n 

<e « w^./j 

- " "^™,0,o(^™(TO(l-e).t)) 

/oo 
d^K{x + iy, t,^ + iT])w{^ + ir))h{^ + irifd^\\L\dx) 
-oo 

(a-i)^ const.e~"* II 1 , ,|2 
<e « 3 w^./i 

Proof. Consider the eigenfunctions equation for the ordinary linear operator Au = u" + p{x) : 

u"+p{x)u = Xu (21) 

Look for a solution exhibiting exponential dichotomy at a; — > ±oo : 

ip+{x, A) = e-T+(^)^a+(x, A), x ^ oo ^^^^ 
(p-{x,X) = e'^-^^^^a-{x,X), x ^ -oo 
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Where, 



7±(A) = a/A-p± 



p± = lim p{x) = (a±) 

X — *±oo 4 

p = max(p+,p_) 



Remark that for | arg(A — p)\ < tt, 



Re(7±(A)) > 



(23) 



(24) 



And consequently, the investigated ip± decrease exponentially as x goes to ±oo. The functions a+ 
and a_ are analytic in x and A and satisfy the following asymptotics : 



a±{x + iy,X) = l + 



1 



, X — > ±oo 



For instance, a+ is a solution as a function of x of the differential equation : 

a'l - + (p- P)a+ = 

which can be reformulated as an integral equation : 

' 1 _ e27+(a:-s) ~ 



a+{x, A) = 1 + 



p 



27+ 



p 



a+(s, \)ds 



This equation may be readily complexified in variable z = x + iy : 



a+{z, A) = 1 + 



P 



27+ 



P 



a+(s + iy, X)ds 



(25) 



(26) 



(27) 



(28) 



On each real line y = const, the function p — p is absolutely integrable and \p{s + iy) — ■p\ds is 
uniformly bounded wrt y as a; — > 00. The integral equation may therefore be solved, which yields 

a a+(2:. A) = 1 + i^\J^\j {a,Yid a a_ satisfying the same estimate) . 

Denote F{x, ^, A) the Green function of the operator A — ^, it satisfies the following estimates : 



\\F{x + iy,i + iri,\)\\L^(d,) < ^e^c.+\i+ixmy-r,) 



|A|3/4 



3_g(C2 + |7+(A)|)(!/-r,) 



(29) 



\\d^F{x + iy,^ + iri,\)\\L'2(dx) < 
Classical spectral theory provides an explicit expression of F{x, ^, A) : 



F{x,i,X) 



1 



W{X) 

1 

W{X) 



ip-{x)ip+{^), if^>x 



(30) 
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Where W(A) = — (7+(A) +7_(A)) +0(1) is the wronskian of the eigenfunctions equation (cf [T5]). 
It is independant of x and ^. Estimates of 11-^(0; + ^ + ir], X\\l2(^^^^ are now requested. 



OO 

■oo 



W{X) 



(31) 



We now estimate the different term as in Sattinger work, using the fact that ||e ^^'^^^\\L'^{dx) 
I^|\y4 , which saves a factor pqiyr compared to L^-norm. For example, consider x > 0, 



|(^_(^ + z77)</.+ (a; + z2;)pd^< e-2M7+(A)):.+2im(7+(A))y f r \^^(^^ + iTj^X)\^d^ 



(32) 



Where, 
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1 / , 1 N \ ^ const 



(33) 



After composition with the square-root, we get the expected power in A {^r[x) estimated by 
Others terms are estimated using the same method. To estimate the norm of derivative, it 



/|A|' 

is sufficient to remark that derivation wrt x of F(x,^, A) gives rise to an extra VA factor. Finally, 
we need to prove that the kernel of the evolution operator e'"^ satisfies the expected estimates and 
that its Ti.^ norm wrt ^ + iry is uniformily bounded in a; + iy. The kernel K{x, t, ^) is represented 
by the following contour integral : 

K{x + iy, t,^ + iT])^ J- I e*^F(x + iy,£,+ iv, X)dX (34) 

Where C is any contour avoiding the essential spectrum of A. Sattinger used the contour A = 
—u! + pe^^^ , parametrized by < p < oo. The two parameter < uj < p and n/2 < 6 < n are 
chosen such that the estimates performed on F{x, t,£^) hold. The same contour will be used. Using 
the estimates on F previously established, we get : 

\\Kix + ty,t,^ + irj)\\L^,^^< const J^l0^e^VWiy-v)\dX\ (35) 

In order to estimate the right handside, rewrite it as an integral on the parameter p = |A + c^l. 

Jl\e^VW{y-v)ux\ < const / . -dp 

|A|3/4 ' '- Jo (pCOsJ-t^)2+p2sin2<5)3/8 

roo gpt cos 5±-ypsin |5| (y— r;) ^ ' 

" Vo C^3/4 + ^3/4|sin^|3/4 
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With f < ^ < TT, cos S <0 and the convergence of the integral is assured. Remark that : 

pt cos S ± ^sm\6\{y - T]) = cos{6) ( Vpi ± 5^ i^^l!?!^ - tan cos 

V COSO 2y/t J 4t ^^^^ 

= -I cos(<5)| (Vjrt± !^ ' + sin \5\^y - 



cos (5 2v^ y ' ' 4t 



Wc conclude by using Young inequality since the kernel K{x + iy, t, (, + iri) is square-integrable with 
respect to x and ^, and the L^(d^)-norm is uniformily bounded wrt x. We use that ^^'ip — 
The estimates are obtained with 7 = 3/4 for if and 7 = 1/4 for its derivative. □ 



4.2 Solving the integral equation 

Proposition 4. Let (X, ||.||) he a Banach Space and B : X ^ X a map such that B{0) = 
and sup < a < 1. Let a G X such that \\a\\ < a. The sequence defined by the 

following iterative scheme : 

hn+i =a + B{hn) 
ho = a 

remains in a ball of center and radius j^oi and has a limit h in X solving the equation h = 
a + B{h). 

Proof. Assume that ||/i„|| < jz^ci- The mean value inequality yields : 

||/in+i|| < a + a\\hn\\ < a + —^—a < — a 

1 — a 1 — a 

The sequence (/i„) is therefore bounded. Using the mean value inequality again, we get : 

— hn\\ < Cr||/l„ — 

Picard's contraction theorem is enough to conclude. □ 

Lemma. The evaluation map ^ b(^c) ^ C is continuous : 



const.{c — \y\) 2 
•j{iy)\i~^°' 



1^(^)1 < , Mnl „ ,(A.) (39) 



Proof. Represent w{z)u{z)^ as an integral using Cauchy's formula over the contour 
C = {^ + ic:^G [00, -00]} U - ic : ^ e [-00, 00]} 

that is to say, 

2 If ^(C)»(C)^rfC .... 
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It provides the estimate : 

1 1 1 

\u{z)\^ < ^^^^ / \w{^ + ic)\.\u{^ + ic)\^d^ (41) 
Vttc- \y\ \w{iy)\ 

which is the expected estimate. □ 
Proposition 5. The nonlinear operator B : Bs,t ^e,T, 



ft roo QQ 

B{h){z,t)= dr —{x + iy,t-T,^ + im{y,T)) 

JQ J~-oo OZ 

R (/(^ + im{y, r)), h{£, + im{y, r))) 



(42) 



with R{f, h) = $(/ + h) - $(/) - ^'{f).h, satisfies the estimate : 

sup \\\DB{h)\\\<Cap{a) (43) 

l[/ilU,T<2a 

With C a constant and p a polynomial of degree < n — 2, both being independant of e and T . 
Proof. Differentiating with respect to h yields : 



OR 
'dh 

where , 



DB{h){g){z,t) = dr J — {x + iy,^ + im{y,T)),t - t) 
-gj^ (/(? + «"z(y, r)), h{^ + im{y, r), r) j g{^ + im{y, r), r)d^ 



(44) 



/c=0 ^ ^ 



For < fc < n — 2, define the following operators : 



f dT 

A'l{h,g){z,t) =- / — - / w{^ + im{y,T))d^K{{x + iy),{^ + im{y,T)),t-T) 
Jo w[z) J_^ 

$('=+2) {f{z))h''+\C + im{y, T),T)g{^ + im{y, T),T)d^ 

1 (It f'^ 

A^{h,g){z,t) = -- — - / + im{y,T))K (x + iy,^ + im{y,T),t - t) 

2 Jo w[z) J_^ 

$('=+2) (/(z))$'(/(z))/i*=+i (e + imiy, r), T)g{C + im.{y, r), r)^^ 

The expected estimates will be proved separately for each and each A^. Proceed first with the 
Aj's : 

/OO 
+ im{y, T))d:cK{x + iy,^ + im{y, T),t- t)^^^+'^^ {f{z))h''hg{^ + im{y, r), T)d^ 

< e^V4eonst|$'=+2(/»)|.|«,(iy)|-(l+«'')'=(i - r)-3/4(yo(l - e) - W'-^ I|/^IIh5,(a„(,„(,_.,.,) 



OO 



X / \w{^ + im{y,T))\.\h{^ + im{y,T),T)\.\g{£,+im{y,T),T)\d^ = h 



OO 
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Where the lemma about the continuity of the evaluation map in H^(ATO(yu(i_e)^^)) and the estimate 
(proposition 2) concerning the kernel dzK have been used. Moreover, the inequality $('=+2) { f{z)) < 
constix;(iy)~''"~'^~^^ hold. The Cauchy-Schwarz inequality leads to : 

h <e*^'/4^onst|«;(iy)|-*("-'=-2)-(i+-^)fe-2a„(i _ r)-^/\y^{l - e) - \y\f-^+''- 



Hence, 



\\\A\{h, .)\\\c(H^) < e^'/*const|«;(iy)r^("-'=-2)-(5+«»)fe-«»-5 



3 fc + 2 

{t - t)4T*(''-i) 



Where, 



And finally, 



/■* dr 
•^0 (t_r)4T'I^^ 



const* = constt'^^'^ 



sup < conste * \iu{iy)\ 

.|U,t<2a 

b„-l , 6, 



L(„_fc_2)-(i+a„)fe-a„-i) 



X (yo - \y\)''^+^e-^H^(^a''+^ 
< consta*=+^e-('"-'^)* 

The last inequality is obtained through the following considerations : 

• The function y 1-^ \'w{iy)\ is bounded on [0, yo] and — ^{n — k — 2) — + an)k — fln — 5 > 0, 



with a„ = - - ^ 



2' 



• The function y^yQ-\y\\s bounded and fc^^^ + > 0, with 6„ = 



• For q > and < tl) < w, t f^e is bounded and q 



n-k-2 
4(n-l) 



> 0. 



Proceed now with A2 

/oo 
wi^ + im{y, t))K{x + iy,^ + im{y, T),t- t)^^^+^^ {f{z))^'{~f{z))h^hg{^ + im{y, r), r)rfC 
-00 

<e^^'^const\^(>'+^\f{z)).<^\f{zmw{iy)\-^^^^^ 

/OO 
\w{(, + im{y, r))|.|/i(C + im{y, r), r)|.|5(^ + im{y, r), T)|d^ = /2 
-00 

Using an analoguous argument as previously, we get : 

h <e^'/4const|«;(it/)|-^(2n-fe-3)-(i+an)fe-2a„(^ _ ^)-l/4(y^ _ \y\f'-^+br. 
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And consequently, 

Jo {t - T)iT^(^^ 

Where, 

dr 1 3n-k-5 

constr U + 47^^-1 = const.i 



1 + 2 

{t - t)4T*<''-i) 



It eventually yields the estimate : 

sup |||4''='(/i,.)|||e,t < const.e^|«;(iy)|-^(2n-fe-3)-(i+«n)fc-«n-ie-c.t^^^„fe+i 

ll'»IU,t<2a 

< const.a'=+ie-('^-^)* 

The differential DB{h){g) being a linear superposition of the monomial operators A\{h,g) and 
^2(^)5), wc got the required estimate. □ 

Definition. Introduce the weighted energy Enit) of a function H{.,t) : 

/oo 
w{0^\H{x,t)fd^ 
-00 

The function H{.,t) is said to have finite weigthed energy if Enit) < 00. 

Remark : A necessary condition for Eh{0) to be finite is that ho{x)dx = 0. 

Proposition 6. The Cauchy problem, for perturbation of Burger's viscous traveling waves has a 
solution h{x,t) which has a holomorphic continuation belonging to B^^t for all e €]0, 1] and t > 
provided that Eh{0) and Eh{0), the weighted energy of the perturbation and of its potential are 
sufficiently small 

Proof. The following integral equation has to be complexified : 

h{x,t)=d,Fo{x,t)- f dr r G{x,t-T;^)R{f{0,H^,r))di (45) 



Assuming that for each r > 0, h{x,T) has a holomorphic continuation in the strip ^m{yQ(i-e),T)- 
The integration contour in the variable ^ can be shifted to the contour ^+im{y, t). The investigated 
h{z, t) solves the resulting integral equation : 

h = d,Fo + B{h) (46) 

where B has already been introduced in Proposition 2. An estimate of the linear term OzFq is now 
required in order to get a control of its norm by some smallness condition on the initial data (here 
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the initial weighted energies Eh{0) and Eh{0))- 

d,Fo{z,t) ^ r ^ (d.K{z,t,0-l:'^'if{z))f{z)K{z,t,o) Hoi^d^ 
J-oo w{z) V 2 J 



OO / -j^ 



<f'(/(0)^(e, ^)u'(0^o(0 + K{z, t, OwiOhoiO- 



Wiz) J-oo V2 

The wave /(^) being bounded wrt ^ € M, using Young inequahty together with the inequahty : 

||if(z,i,e)||LHd?) < const e*^'/4 

we get the following estimate of dzFo : 

V \w{iyo{l - £))|^^ / kWo(l - e))\y^ 

Proposition 1 provides a solution if 

Cap{a) < 1 (48) 

□ 

Proof of theorem 1 : 

One can use the following a priori estimates for the weighted energy and conclude that initial data 
may actually be considered as always sufficiently small 

Proposition 7. If Eh{0) is finite, then Enit) and Eh{t) are finite for all t > and are bounded 
from above by Ciexp {—C2t). 

This statement is not well-known, but it is a variation of the result of Iljin-Oleinik [7], its proof 
can be found with all details in [16| . Now, solutions of our main theorem may be constructed 
through the procedure : '"wait"' until (|48|) is fuUfillcd denote this date T* , and apply Proposition 
3 with h{x,T*) as initial data. 
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